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INTRODUCTION 


Several random processes have been proposed as mathematical models 
for atmospheric turbulence in aeronautical applications- The amplitude 
modulated random process or the Press model (references 1 and 2} has been 
■widely applied to the analysis and measurement of aircraft response to 
atmospheric turbulence (references 3 to 7> for example). The amplitude 
modulated process accounts for the local variations of the turbulence in 
combination with a random modulation of the amplitude. In reference 8 this 
process was extended by the addition of an independent mean value process, 
based upon a reinterpretation of the random process described in references 
9 to 12. The resulting amplitude-modulated-plus-mean (AMPM) random process 
allows more versatility in modeling the measured properties of atmospheric 
turbulence. The basic mathematical properties of the AMPM process, including 
the first and second order properties and the characteristic function of 
general order, have been developed in reference 8. 

The correspondence between the properties of the AMPM process and 
the associated properties of both atmospheric turbulence and aircraft 
response are examined on a qualitative basis in the present study, which 
consists of three parts. The first part is an interpretation of the _AMPM 
random process in terms of the physical structure of atmospheric motions. 

The second part is an examination of the considerations involved in 
application of the AMPM process as a model for atmospheric motions. The 
third part is an evaluation of the effect of the random mean value 
variation upon aircraft response. A procedure for the development of 
aircraft strength design criteria is outlined. The mathematical properties 
of the AMPM process are reviewed in Appendix A. 
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standard deviation of either the subscripted process or 
(without subscript) the R process 
filter constant of aircraft plunging motion 
standard deviation of S process 
amplitude intensity parameter 
standard deviation of M process 
mean value intensity parameter 
mean aerodynamic chord length of wing 
ensemble average 
error function (reference 13) 
complementary error function = 1 - erf ( ) 
scale of turbulence of subscripted process 
fourth order flatness factor or kurtosis ( fourth moment 
~ square of second moment) 
mean value random process 

expected frequency of positive slope crossings of indicated 
level 

expected frequency of positive slope zero crossings of R 
process 

probability parameter , 

local random process 

amplitude random process 

integral scale of subscripted process 

time 

gust intensity parameter 
aircraft speed 
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amplitude-mo dul ated-plus-mean ( MPM ) random prucess 
aircraft plunging velocity 

vertical velocity component of atmospheric motion 

deterministic mean value of w 

intensity random variable 

amplitude modulated random process 

ratio of standard deviations of Z and M processes 

ratio of N of subscripted process to N 

o or 

aircraft mass parameter 

autocorrelation function of subscripted process, correlation 
coefficient of subscripted variables 
standard deviation of subscripted process 
time difference variable 

power spectral density function of subscripted process 
normalized power spectral density function of subscripted 
process, 0 <_ (P <_ 00 

autocovariance function of subscripted process 
frequency, Fourier transformation variable of T 
natural frequency of dynamic system 

breakpoint 
input random process 
output random process 
quasi- steady 



PHYSICAL INTERPRETATION OF THE AMPM PROCESS 


The AMPM process is formed "by the sum of an amplitude modulated and a 
mean value random process. The amplitude modulated process is formed hy the 
product of a local R component process and a slowly varying amplitude S 
component process which modulates the local process. 

z(t) = r(t) s{t) (1; Al)* 

The amplitude modulated process accounts for some of the Gaussian and 
non-Gaussian aspects of atmospheric turbulence measurements. For a short 
sample of the amplitude modulated process, the slowly varying amplitude 
component is approximately constant; measurements of the process have 
approximately Gaussian properties. For a long sample, the amplitude 
component varies as a random function of time; measurements of the process 
have strongly non-Gaussian properties. The defining relations for the 
AMPM (amplitude-modulated-plus-mean) or total process are 

w(t) = z(t) + m(t) (2; A2) 

w(t) = r(t) s(t) + m(t) (3; A3) 

*The dual equation numbers identify equations which are discussed in 
Appendix A. Equation (l; Al) is both equation (l) of the text and equation 
(Al) of Appendix A. 
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The three component processes are identified as the local R, the amplitude 
S, and the mean value M. * Following the development of references 9 to 12 
the three component processes are specified to he independent, stationary, 
and Gaussian with zero mean values. 

A physical interpretation of the AMPM random process, which considers 
both recorded atmospheric turbulence data and general turbulence properties , 
is presented in this section. Two aspects of the random process are of 
primary interest. The first is the addition of the independent mean value 
component to the amplitude modulated process.. The second is the interpre- 
tation of the mean value component as slowly varying relative to the 
amplitude modulated or, more specifically, relative to the local R componenl 
process. This aspect is a reinterpretation of the original development of 
references 9 to 12, which assigned the same integral scale value to the 
three component processes. 

The qualitative effect of slowly varying changes in both the amplitude 
and the mean value can be evident in measured atmospheric turbulence data. 
Figure 1 shows an example of measurements of the three velocity components 
of atmospheric motion. The data show two periods of moderate turbulence 
superimposed on a background of light turbulence. These two periods show 
significant increases in the amplitude of the local variations. Slow, 
random variations of the mean value are also apparent. The vertical velocity 
component shows a dominant amplitude modulation with relatively little 

*The term mean value process is used for m(t ) , since it appears as a 
random variation in the mean value of the local R process. Care 
must be used to distinguish between this mean value process and the mean 
value of any measurement of the AMPM process. 
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variation of the mean value. The two horizontal components, particularly 
the lateral, show large, slow variations in the mean value. 

The set of data in figure 1 is a case of clear air turbulence related 
to high-altitude wind shear. The data were recorded above the Death Valley 
area of California at an altitude of 13 km with an airspeed of 183 m/s. 

The average wind direction was almost perpendicular to the flight path. 

The data were obtained as part of the MAT (Measurement of Atmospheric 
Turbulence) program at HASA-Langley Research Center (reference 1^+ ) . 

Special care was fcaken in the measurement of the long wavelength part of 
the atmospheric motion, particularly in removing the low frequency aspects 
of the aircraft motion (reference 15)- Thus the slow variation of the mean 
value shown in figure 1 is atmospheric motion; the residual aircraft motion 
in the data is negligible. 

The properties of the AMPM process can be interpreted by examining 
the associatevd power spectral density function, which is the sum of the 
functions of the independent amplitude modulated and mean value processes. 

4> (u) = $ (w) + $ (oj) (U; A23) 

w z m 

The power spectral density function consists of two terms with significantly 
different frequency dependence. The composition of the resulting power 
spectral density function is shown in schematic form in figure 2. The 
term corresponding to the mean value process has a predominantly low 
frequency content, since this component is slowly varying with respect to 
the local K component process. The term corresponding to the amplitude 
modulated process has a predominantly high frequency content, which is 
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associated with the local R component process. (The effect of the 
amplitude process upon the spectral properties of this term is relatively 
small as will he shown subsequently). The higher frequency term, 
corresponding to the amplitude modulated process, is currently accounted 
for in aeronautical applications. The lower frequency term is an 
additional effect which is included in the AMPM process. 

The structure of the AMPM process suggests a possible correspondence 
witn the basic composition of atmospheric motion, which consists of several 
different types of motion (references 1 6 and 17, for example). The primary 
classification divides atmospheric flow into two basic kinds of motion: 
internal gravity waves (winds or drafts) and turbulence. Based upon the 
associated ranges of frequency values, the internal gravity waves would 
correspond to the lower frequency mean value process and the turbulence 
■would correspond to the higher frequency amplitude modulated process. 

The resulting separation of the associated power spectral functions into 
two terms in the manner of figure 2 appears in studies of measured 
atmospheric motions (references 18, 19 and 20, for example). A specific 
example of the separation of the spectral function into two distinct terms 
is given in reference 21, which includes an examination of the transfer 
of energy between the terms. 

The possible correspondence between the AMPM process and the composition 
of atmospheric motions must be thoroughly examined before it is established 
on a firm basis. Such an examination must consider the basic mathematical 
properties of the random process and the corresponding physical properties 
of atmospheric motion. In particular, the random process specifies two 
important properties. First, the AMPM process specifies that the two 
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kinds of atmospheric motion have significantly different probabilistic 
structures, the lower frequency term being Gaussian and the higher frequency 
term being formed by the product of two independent Gaussian processes. 

This difference is reflected in all associated properties, for example, the 
probability density function and the associated flatness factors. Studies 
such as references 18 to 21, which show two terms in the spectral functions, 
do not examine the probabilistic structure associated with these terms. 
Second, the AMPM process Specifies that the two terms are statistically 
independent. Since the two terms consequently can not interact, the AMPM 
process can not account for the interaction due to nonlinear coupling 
between the two kinds of motion, which is an essential physical property 
(reference 1 6 , for example). One possible means of removing the 
independence property is the introduction of correlation between the 
amplitude and mean value component processes. This modification and the 
properties of the resulting AMPM process are discussed in appendix B, 
where it is shown that the effects of the correlation upon the probabilistic 
structure are small if the amplitude and mean value processes are weakly 
correlated. 

An associated question is the determination of the values of the 
physical parameters associated with the random mean value variation in 
atmospheric motions. Although the answer to this question requires an 
extensive examination of atmospheric data, a preliminary estimate is made 
to determine representative values for use in the subsequent discussion. 

The first quantity of interest is the a parameter, which is the ratio 
of the standard deviations of the amplitude modulated and the mean value 
processes. 
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(5; Al4) 


In reference 9 a value of slightly less than one is suggested for the a 
parameter of the vertical velocity of atmospheric motion, based upon 
examination of measured data. The value of one is used as a representative 
value of the a parameter in the subsequent discussion.* The second 
quantity of interest is the ratio of the integral scale values of the mean 
value and the amplitude modulated processes. Based upon published 
spectral data which show the existence of both processes (references 18 
to 2l), the value of ten is used in the subsequent discussion as a 
representative value of this ratio for the vertical velocity component of 
atmospheric motion. 


*If the a parameter is equal to one, then b = c since A is given 
a unit value for turbulence processes. The condition of equal Standard 
deviations (b = c) does not imply that the values of the amplitude and 
the mean value processes are equal within any short segment of turbulence, 
since those values are only member functions of their total ensembles. 



APPLICATION OF THE AMPM PROCESS AS AN ATMOSPHERIC TURBULENCE MODEL 

The considerations involved in modelling the properties of atmospheric 
motions with the AMPM process are examined in this section. The primary 
topic is the effect of the slowly varying mean value component , since the 
presence of this component is the difference between the amplitude modulated 
and the AMPM process. The discussion begins with an examination of the 
correlation and spectral properties of the AMPM process. The relationships 
between these dynamic properties and the probabilistic structures of both 
the AMPM process and atmospheric motions are then examined. The effect of 
the mean value -Variation upon the exceedance expression of the AMPM process 
is also examined. 

Correlation and Spectral Properties 

The mean value variation has a strong effect upon the autocorrelation 
function of the AMPM process, particularly for large values of the time 
difference. There are corresponding effects upon the power spectral 
density function, particularly for small values of the frequency. In 
the present understanding of atmospheric turbulence properties there 
are open questions on the spectral properties at low frequency values 
and consequently on the value of the integral scale (reference lU) . 

The (normalized) autocorrelation function of the AMPM process is 
determined from the corresponding functions of the three component _ 

processes > 

P (t) = -~r P z (t) + " 2 -- — P w (t) (6; A3 la) 

or + 1 a + 1 
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where p a (x) = p r (x) P s (t). 

The corresponding relation for the normalized power spectral density 
function is the Fourier transformation of equation (6). 

2 

<(> (m) = * (w) + * (to) (7 ; A3Ib) 

w a 2 + 1 2 a 2 + 1 m 


The preceding relations show the relative contributions of the amplitude 
modulated and the mean value processes to the functions of the AMPM process. 
The relative contributions are determined by the a parameter. For small 
values of the a parameter, the functions of the mean value component are 
dominant; for large values of the a parameter, the functions of the 
amplitude modulated process are dominant. The autocorrelation and spectral 
functions of the AMPM process must be distinguished from those of the 
three component processes. Most discussions in the aeronautical literature 
and most applications to aircraft design (references 3 to 6, for example) 
consider only the functions of the local R component process. 

One measure of the dynamic properties of a random process is the 
integral scale, which is the integral of the autocorrelation function 
over the range of positive time values. 


T 
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w . . 

- / P (t) dT 
o K w 



The integral scale is related to the zero frequency value of the normalized 
power spectral density function, the relation following from equation (8) 
and from the definition of the Fourier transformation. 
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T = $ (0) 

w 2 v' 


(9) 


These relations for the integral scale follow the notation used in the 
atmospheric turbulence and the aeronautical literatures (reference 5, 
for example). However, the autocorrelation is usually defined as a 
function of the spatial dimension; the corresponding integral scale is 
referred to as the scale of turbulence. Also, the present notation uses 
W to denote any of the vector components of the atmospheric velocity, 
not specifically the vertical component. 

The integral scale value of the AMPM process is determined by those of 
the amplitude modulated and mean value processes. The relationship is 
obtained by combining the definition of the integral scale, equation (8), 
and the relation between the autocorrelation functions, equation (6). 
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In the concept of the AMPM process, the mean value process is slowly 
varying relative to the local R component and also relative to the 
amplitude modulated process. The value of the integral scale of the 
mean value process is accordingly much greater than that of the amplitude 
modulated process. Thus by equation (10) the presence of the mean value 
variation strongly influences the- value of the integral scale of the 
AMPM process. This can be true even if the contribution of the mean 
value component to the standard deviation of the total process is fairly 
small, that is, if the value of the a parameter is large. 



The amplitude component process also influences the integral scale 
of the AMPM process, although to a much smaller extent. This effect 
can he seen by a specific example. The functions of both the local R 
and the amplitude component processes are assumed to have the form for 
the longitudinal velocity component in the Dryden model for turbulence 
(reference 5 ). Equivalently, the Gaussian R and S components are 
assumed to be Markov processes. Their autocorrelation functions are 
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The autocorrelation function of the amplitude modulated process is the 
product of the autocorrelation functions of the local and amplitude 
component processes by equation (6). 

-M/t z , s 

P z (t) = e (12) 

•where + ~- 

z r s 

Under the modulation concept the integral scale of the amplitude process 
is much larger than that of the local R process. In this case, the 
amplitude modulation has the effect of slightly reducing the integral 
scale value of the original R process. For example, if the integral 
scale of the amplitude process is ten times that of the local R 
process, then the amplitude modulation reduces the integral scale of 
the Original process by about nine percent. The amplitude component 
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accordingly has a minor effect on the correlation and spectral functions of 
the amplitude modulated process. Conversely, it will he difficult to 
separate the autocorrelation and spectral functions of the local R and 
the amplitude components from the measured functions of the amplitude 
modulated process. The preceding example is a special case where the 
separation is impossible. 

Effect of Dynamic Properties Upon Prohatilistic Structure 

The mean value variation has significant effects upon the probabilistic 
structure of the AMPM process. These effects are related to the dynamic 
properties of the process, since the AMPM process is the sum of two 
independent processes which have significant differences in both their 
dynamic properties and probabilistic structures. 

As an example of the effects of the mean value variation, the 
differences between the AMPM random process and its first derivative are 
considered. The defining relation for the AMPM process is 

w(t) = r(t) s (t ) + m(t) (3) 

The quasi-steady derivative of the process is considered.* 

:: w_(t) — r(t) s(t) (13; A50) 

*The quasi-steady form of the AMPM process is used in the discussion, 
that is, the dynamic properties, such as the derivatives, of the slowly 
varying amplitude and mean value components are omitted. The dynamic 
properties of the local R component are retained. 
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For the quasi-steady derivative the value of the a parameter is infinite. 
Consequently the AMPM process and its first derivative generally have 
different relative contributions from the amplitude modulated and the mean 
value processes 3 that is, they have different values of the a parameter. 
The difference is shown by the values of the associated fourth order 
flatness factors. 


Mi. (w) 


U p 

3 (3a + 2a + 1) 
(a 2 + l) 2 


(14; Al6) 


M 4 (w ) = 9 


(15) 


The flatness factor of the derivative follows from the general form of 
equation (l4) with an infinite value of the a parameter. The value of 
nine also follows from the form of the derivative, equation (13), which is 
the product of two independent Gaussian processes. 

The amplitude modulated process and its quasi-steady first derivative 
have similar properties. 

z(t) = r(t) s(t) (l) 

z (t) = r(t) s(t) (l6) 

The amplitude modulated process and its quasi-steady derivative have 
the same basic form; the product of a Gaussian local process and the 
same amplitude process. Consequently their probability density functions 
have the same functional form and their fourth order flatness factors 
have the same value. 
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M 4 (z) = M u (£ ) = 3M 4 (s) ( 17 ) 

The flatness factor has a value of nine if the amplitude process is Gaussian. 
Other distributions have been suggested for the amplitude process, enabling 
the resulting amplitude modulated process to account for a range of values 
for the flatness factor. However, in all cases the resulting amplitude 
modulated process and its quasi-steady derivative must have the same value 
for the flatness factor. 

vonsequently one effect of the slow, random mean value variation is 
that it allows the AMPM random process and its derivative to have different 
probabilistic structures. This property cannot be accounted for by 
amplitude modulation alone. There are several studies which show that this 
property appears in measurements of atmospheric motions. These studies 
present a strong argument that the inclusion of the mean value variation 
gives a better representation of the properties of atmospheric motions. 

One study is reference 22, which examined turbulence related to storms. The 
probability distribution functions of the atmospheric velocity components 
and their gradients were examined and presented in graphical form. Comparison 
of the measured probability distributions shows that the functions of 
the velocity components and their gradients have significantly different 
functional forms. The gradients show a much greater deviation from the 
Gaussian form, a property which is consistent with the form of the flatness 
factors of the AMPM process, equations (14) and (15). 

A more definitive study is that of Chen, reference 23. Atmospheric 
turbulence data were obtained from several sources: low altitude. 



severe storm, and high, altitude clear air turbulence. The properties 
of the three velocity components and their gradients were examined, 
including computation of the fourth order flatness factors. The results 
show a consistent trend of higher values of the flatness factor for the 
gradients than for the velocity itself. The flatness factors for the 
velocity generally had values somewhat greater than three, with a few 
values being slightly less than or equal to three. The flatness factors 
for the gradients had values between three and seven. These values are 
less than the value of nine for the derivative of the quasi-steady process 
However values below nine are possible if transition effects* are present 
or if the quasi-steady approximation is not completely satisfied, 
particularly since the results are based upon the gradients and not the 
derivatives of the velocity components. Thus the results of reference 23 
are generally consistent with the properties of the AMPM process. 

The study of references 2k and 25 presents data for atmospheric 
motions which were measured by tower-based instrumentation. The 
probability density functions and the associated flatness factors were 
determined. The data show higher values of the flatness factor for the 
velocity gradients than for the velocity itself. Another set of results 
was obtained in an examination of the effects of high-pass filtering 
upon the measured data. The filtering consisted of removing a running 
mean value from the data. The filtering calculation considered several 
time periods , the shorter time periods removing more of the lower 

^Transition effects are related to the development of the slowly varying 
amplitude and mean value components from their initial values in any 
measurements of the AMPM process (reference 8). 
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frequency content of the velocity components. Using the concept of the 
AMPM process, this operation represents an approximate hut direct removal 
of the random mean value variation. The resulting flatness factors show 
a definite trend toward higher values as the filtering operation removes 
more of the low frequency content. These results indicate that the 
inclusion of a slow mean value variation gives a letter representation of 
the measured atmospheric motions. The highest values of the flatness 
factor were generally less than the value of nine for the quasi-steady 
derivative, which represents the complete removal of the mean value 
variation. The data also show some values of the flatness factor less than 
three and some relatively large values for the skewness (normalized third 
moment). These results can not he modelled by the AMPM process in its 
fully developed form, hut could represent the transitional form of the 
process (reference 8). 

The difference between the prohahilistic structures of the velocity 
and the velocity gradient of atmospheric motions has been examined in 
other studies where possible explanations of this property have been 
suggested. For example, reference 2 6 uses this property to recommend that 
the gradient rather than the turbulence velocity is the more important 
experimental quantity. The turbulence velocity is considered to be a 
Gaussian process, thus omitting the amplitude modulation effect. The 
analysis of reference 26 shows that the Gaussian property is lost in taking 
the gradient of the process. However the gradient, which is the difference 
between the values of the process at two time values, is a linear 
combination of Gaussian random variables and therefore must itself be 
Gaussian (reference 27). Consequently the development of reference 26 


can. not explain the difference between the probabilistic structures of 
the velocity and the velocity gradient of atmospheric motion. 

In summary, the probabilistic structure of the AMPM random process 
possesses some mathematical properties that are not present with 
amplitude modulation alone. In particular, the probabilistic structure 
of the AMPM process can vary as different ranges of frequency values are 
emphasized. Measurements of atmospheric motions show corresponding 
properties, for example, the velocity gradients show more deviation from 
Gaussian properties than the velocity itself. Atmospheric motions thus 
shov significant properties which can be accounted for by the AMPM process. 
It is noted that this conclusion is in disagreement with that of 
references 11 and 12, which are the original development of the basic 
random process. However the development of those references did not 
include the differences between the dynamic properties of the three 
component processes, which are essential for the existence of the 
indicated dynamic properties of the AMPM random process. 

Exceedance Expression 

The mean value variation has significant effects on the exceedance 
expression of the AMPM process. An analytical relation for the exceedance 
expression has been developed by using the quasi-steady approximation. 

The relation is expressed as an exceedance ratio: the ratio of the 

expected frequency of crossings of a given level of the AMPM process to 
the expected frequency of crossings of the zero level for the R 


component process. 
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i«M = i e 1/2 “ 2 . {e- w/ac [l + erf (— 1-)] 

or 2 JSe 

+ e w/,ac erfc(-^— + )} (l8; A^k) 

/2c /2a 

where w >_ 0 

N;(-w) = N(w) 

The exceedance ratio is plotted in figure 3 as a function of the ratio of 
the level of the AMPM process to its standard deviation, and for several 
values of the a parameter.; For large values of the a parameter the 
exceedance ratio approaches the exponential form of the amplitude modulated 
process, equation (A 5 T)* For small values of a the exceedance ratio 
approaches a Gaussian form related to the mean value component process, 
equation (A58). The intermediate cases show a combination of those two 
functional forms. At the zero level the exceedance ratio is not equal 
to one since the exceedances are ratioed to the expected number of zero 
crossings of the R component process and not that of the AMPM process. 

The exceedance ratio can he strongly influenced hy the quasi-steady 
approximation at low levels of the process and for small a values, since 
the dynamic properties of the mean value component are absent. The exact 
form of the exceedance expression is examined in appendix C. 

The contribution of the mean value variation to the exceedance ratio 
of the AMPM process is not clearly shown in figure 3 , since the mean value 
process affects the non-dimensional level of the AMPM process through its 
standard deviation, which depends upon the a parameter. 
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a 2 = c 2 (a 2 + 1) 
w 


(19; A15) 


The contribution of the mean value variation is more clearly shown in 
figure where the exceedance ratio is plotted as a function of the ratio 
of the level of the AMPM process to the standard deviation (Ab) of the 
amplitude modulated process, which is independent of the a parameter. 
Figure it thus shows the effect of the mean value variation relative to the 
amplitude modulated process, which corresponds to an infinite value of. the 
a parameter. Figure it shows that the mean value variation can either 
increase or decrease the exceedances depending upon the values of both the 
process level and the a parameter. The mean value variation can 
significantly increase the exceedances at large values of the process as 
the value of the a parameter approaches zero. In the limit of large 
values of the process, the exceedance ratio approaches the exponential 
form of the amplitude modulated process, but with the exceedances increased 
by a constant factor whose value follows from equation (A59). 


lira. N(w; a) _ l/2a 2 

w“*» F(w; a=°°) e 


( 20 ) 


This factor indicates the upward shift, with decreasing value of a, of 
the limiting exponential form of the exceedance ratio curves of figure k. 

*» The effect of the mean value variation upon the exceedances can be 

i seen in the measurements of the velocity of atmospheric motion in figure 1. 

I,'-- *■' . 

The lateral velocity component shows crossings at high levels which are 

i • - • ' ■•■■■■ 

[' ' ■ " ' • - ... . ■ . • . ‘ .-. •• ,' ■■ ■ ■' •' • • 

related to the combination of the rapid local variations superimposed on 
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the large, slow variations of the mean value. For example, the 
crossings of the -20 m/s level are related to the mean value variation. 
Crossings of this level do not occur in the other two velocity components 
which show similar local variations hut not the large variations in the 
mean, value. Also, the mean value variations cause the crossings at the 
zero level of the lateral velocity component to he much less than those 
of the other two velocity components. The measured data thus show the 
decreased exceedances at low levels and the increased exceedances at high 
levels, due to the mean value variations, which are indicated in figure 4. 
Experimental Aspects 
Intensity Process 

It was estimated previously that the standard deviations of the 
amplitude modulated and mean value processes are approximately equal for the 
vertical velocity component of atmospheric motion. The condition of equal 
standard deviations for these two processes suggests an explanation for 
another aspect of measured atmospheric turbulence data. Attempts have been 
made to determine the distribution of the "intensity" process directly 
from measured data, in order to verify the assumed Gaussian distribution of 
the amplitude process. Reference 28, for example, shows that the measured 
data suggest a Rayleigh rather than a Gaussian distribution for the 
intensity. The measured data were obtained from two studies, references 
29 and 30, which involved direct measurement of atmospheric motions. The 
resulting intensity random variable is 

x - \ S 2 + m 2 . . : (21) 
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vher e E[r^] = 1 

This assumes that the data samples are of appropriate length to allow the 
full development of the local R component without significant development 
of the amplitude and mean value component processes * The resulting intensity, 
x, is a random variable which depends upon the values of both the amplitude 
and the mean value processes. Following the basic formulation, the random 
amplitude and mean value variables are independent and Gaussian with zero 
mean values. If they also have equal standard deviations, the intensity x 
of equation (21) has a Rayleigh distribution (reference 2 r f)- Thus the 
distribution of turbulence intensity indicated by the results of reference 
28 suggests a correspondence with the properties of the AMPM random process. 
Experimental Procedures 

The determination of the spectral properties of the mean value 
variations in atmospheric motions presents a difficult experimental problem. 
There are several problems in the accurate measurement and analysis of 
atmospneric motions in the range of low frequency values associated with 
the mean value process, especially when the data are measured with a 
flight vehicle. References 5 (appendix C), 15 and 31 give discussions of 
the difficulties associated with the removal Of the extraneous effects of 
aircraft motion from the measurements. Another problem is the accurate 
computation of the spectral functions. Reference 32, for example, shows 
that significant spurious effects can be introduced in the low frequency 
range by the numerical procedures used to compute the spectral functions. 

Thus caution must be used in both the measurement and the analysis of 
atmospheric motions in order to determine the spectral properties of the 
mean value variation. 



EFFECT OF THE MEM VALUE VARIATION UPON AIRCRAFT RESPONSE 


The effect of the mean value variation in atmospheric motions upon 
aircraft response is examined in the present section, using estimates of 
the dynamic properties of Loth the mean value process and significant 
aircraft response quantities. Also, a procedure for developing design 
criteria for aircraft structural strength is outlined. 

AMPM Process and Aircraft Response 

The effect of the AMPM process upon aircraft responses can be examined 
by considering a simple mathematical model, in which the aircraft is a 
rigid body in level flight at constant airspeed and is allowed to move in 
plunging motion only. The corresponding equation of motion is 

w + aw = aw (22) 

g 

where w is the plunging velocity of the aircraft, w the plunging 

acceleration and w is the vertical velocity of the atmospheric motion. 

g 

The plunging velocity has the characteristics of a low-pass filter; the 
frequency response function is approximately constant for frequency values 
below the filter constant a and decreases as a. minus-one power of 
frequency above a. The acceleration has the characteristics of a high- 
pass filter; the frequency response function increases proportionally to 
frequency below a and is approximately constant for frequency values 
above a. The filter constant' of the first-order system is related to 
aircraft response parameters. •; 
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where V is the aircraft speed, c the mean aerodynamic chord of the wing, 

and y the dimensionless aircraft mass parameter (reference 33)*. 

6 

The aircraft plunging motion is developed for the response to atmos- 
pheric motion which is modeled by the AMPM process. The spectral properties 
of both the local R and the mean value components are taken to be those 
of the transverse velocity component in the Dryden model of turbulence 
(reference 5). 
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where to , = 1/T r = V/L and L is the scale of turbulence of the R 
rb r r r 

component process. For the first-order differential equation, equation (22), 
and the normalized spectrad. density function for the transverse velocity 
component, equation (24), the variances of the response are (reference 34) 
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The aircraft plunging response to the AMPM process is developed from 
the previous relations by using the quasi-steady approximation for the 


*\i fS =■ 2(W/S)/(gc&k p) , where (W/S) = wing loading, g = gravity constant, 

C = lift curve siope, and p = air density. 

L a 
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amplitude modulated process. The variance of the response to the AMPM 
process is the sum of the variances for the responses to the independent 
amplitude modulated and mean value processes. The a parameters of the 
response are determined by the standard deviations of the response to the 
local Tl and mean value components and by the a parameter of the AMPM 
process, for example. 





(26; A48) 


The a parameters of the plunging response are determined by combining 
the previous relations and introducing the appi'opriat e notations for the 
two component processes.* 
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(27b) 


The a parameters of the response quantities depend upon the three 
frequency constants: that of the plunging response, a, and those of the 

two component processes, to ^ and to^ . 

The preceding analysis is used to estimate the effects of the mean 
value variations in atmospheric motions upon the aircraft plunging motion 


*The dynamic properties of the mean value process have been included in 
the variances of the response, but have been omitted in the exceedance 
expression, equation (l8) . The quasi-steady values of the variances are 
obtained by setting to to zero. 
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for the present class of subsonic transports. The ratio of the frequency 
constants of the plunging response and the local R component process is 



(28) 


Using typical parameters for subsonic transport aircraft, the mass parameter 
value is usually between 50 and 100. Using typical values of the wing chord 
length and the scale of turbulence, the ratio of the frequency constants has 
a minimum value of about one, that is , a = w . (Larger values of the ratio 
of equation (28) result in larger values of a and thus less effect of 
the mean value variation upon the plunging acceleration.) Using the previous 
estimates of the physical properties of atmospheric motions, the a 
parameter of the vertical velocity is equal to one, and the ratio of the 
integral scales of the mean value and the local R components is equal 
to ten. With these values of the parameters, the resulting values of the 
a parameters and the fourth order flatness factors of the plunging 
response are 


M 1+ (w; a - .66) = 3.55 

M 4 (w g ; a = 1.00) = 4.50 (29)* 

M 1( (w; a =. 2.18) = 7.09 


*If the quasi-steady approximation is used in determining the variances 
of the response to the mean value process, that is, w ^ is set equal to 
zero, the results are M^(w; a = . 6l) = 3.44 
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The results show the effects of the dynamic system response upon the relative 
contributions of the lower frequency mean value process and the higher 
frequency amplitude modulated process. For the plunging velocity, which has 
the characteristics of a low-pass filter, the relative contribution of the 
mean value process is increased, resulting in a lower value of both the a 
parameter and the fourth order flatness factor. The opposite effects are 
shown by the plunging acceleration, which has the characteristics of a 
nigh-pass filter. 

The associated effects upon the quasi-steady exceedance ratios for 
the three quantities of equation (29) are shown in figure 5 > with the process 
level ratioed to the standard deviation of the response to the amplitude 
modulated process in the same manner as figure k. For the original AMPM 
process (a parameter equal to 1.00) the effect of the mean value 
variation is small. For the plunging acceleration (a equal to 2.18) the 
effect is almost negligible. For the plunging velocity (ot equal to 0.66) 
the mean value variation significantly increases the exceedance ratio, 
except at the lowest response levels. The results of figure 5 show two 
points. First, with the presence of both the amplitude modulated and the 
mean value processes, the exceedance expression depends upon the dynamic 
properties of the system response quantity and consequently can be different 
for different response quantities. Second, the presence of the mean value 
variation can significantly increase the exceedances, particularly for a 
response quantity which acts as a low-pass filter. Since the critical 
aircraft response quantities are usually either accelerations or incremental 
loads which are closely related to accelerations, the effect of the mean 
value variations in atmospheric motions upon aircraft response are 
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generally not significant. The effect of the mean value variation however 
can he important for response quantities which are related to the velocity 
or displacement. 

The effects of the mean value variation in atmospheric motions upon 
aircraft response are also examined by considering quantities which have 
the characteristics of a narrow hand-pass filter, such as low damped 
oscillatory modes, either whole body or elastic. For simplicity both the 
amplitude modulated (or equivalently the local R) and the mean value 
processes are assumed to have the Dryden spectral function for the 
longitudinal component of atmospheric turbulence (reference 5), for 
example , 
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The power spectral density function has a simple, approximate form: for 

frequency values below Oi^ the function is constant, and for values above 
ui^ the function decreases as the minus-two power of the frequency. This 
approximate form is plotted in figure 6 for both the mean value and the 
amplitude modulated (or the local R) processes. Since the mean value 
component is considered to be slowly varying, the integral scale value 
is significantly greater than that of the local R process, giving the 
opposite relation for the corresponding values of the frequency breakpoint 
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The power spectral density function shows the basic effects of the 
two random processes upon the response of dynamic systems which have the 
characteristics of narrow band-pass filters. The associated frequency 
response function is dominated by a single frequency value, which is the 
natural frequency 03 n of the dynamic system. The relative effect of the 
two random processes upon the variance of the response is directly 
related to the values of the two power spectral density functions at the 
value of the system natural frequency. The effect can be visualized 
from figure 6. The functional properties of the composite power spectral 
density are divided into three regions of frequency values . First , for 
values of the system natural frequency which are below the frequency 
breakpoint of the mean value process , both spectral functions are constant 
The response is dominated by the mean value process; the a parameter of 
the system response is smaller than that of the excitation process. 
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Second, for values of the system natural frequency which are between the 
two frequency breakpoints, the spectral function of the amplitude 
modulated process is constant while the function of the mean value 
process decreases with increasing frequency. Either process can be 
dominant in the system response, depending upon the value of the natural 
frequency. 
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Third, for values of the system natural frequency which are above the 
frequency breakpoint of the amplitude modulated process, both spectral 
functions decrease as a minus-two power of frequency. The response is 
dominated by the amplitude modulated process; the a parameter of the 
system response is larger than that of the excitation process. 
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Consequently the relative contributions of the two processes to the response 
of a narrow band-pass system to the AMPM process depend on the value of the 
natural frequency of the dynamic system relative to the values of the 
frequency breakpoints of the amplitude modulated and the mean value processes. 

The preceding development can be used to estimate the effects of the 
mean value variation in atmospheric motion upon aircraft dynamic systems. The 
frequency breakpoint of the amplitude modulated process, which is 
approximately equal to that of the local R component process, can be 
estimated from the values of the scale of turbulence currently suggested in 
the aeronautical literature. There is an open question on the most appropriate 
value, but estimates vary from about 250 m (800 ft, reference 35) up to 
the value of 762 m (2500 ft) used in aircraft design criteria, references 3 
and k. For the present class of subsonic transports this gives values of 
of between about 0.2 and 1.0 radians per second.* The dominant whole body 


*The frequency breakpoint is given by the relation w 


L 5 


where the constant 


r 

y locates the "knee" of the spectral curve and has a value between 1.0 and 1.5 
for the atmospheric turbulence models used in aeronautical applications. 



32 


modes of aircraft, the Dutch roll (lateral) and the short period (longitudinal) 
modes, have representative values of about one and two radians per second, 
respectively, for subsonic transport aircraft. The values of the primary 
response frequency thus fall mostly in the third region of the spectral function, 
which is dominated by the amplitude modulated process. Consequently, the value 
of the a parameter for most aircraft response quantities is greater than 
that of the original atmospheric velocity components. 

The effects of the mean value variation in atmospheric motions are 
thus fairly small for most aspects of the response of current subsonic 
transport aircraft, for example, response quantities such as accelerations and 
incremental loads which are related to acceleration quantities. The mean value 
variation will be more important for response quantities related to the velocity 
or displacement. The mean value variation may significantly effect the long 
period or phugoid mode associated with speed and altitude perturbations of 
aircraft (reference 36). This mode is usually ignored in the calculation 
of aircraft loads. However, there are cases where the effects of the phugoid 
mode may be significant. The effects of the mean value variation may become 
important for aircraft which are significantly larger than the current class 
of subsonic transports, due to the generally lower values of the frequencies 
of the whole body modes. The mean value variation is more important for 
higher flight speeds due to the corresponding higher values of the frequency 
breakpoint of the amplitude modulated process. A similar conclusion applies 
to flight at low altitudes due to the associated lower values of the scale of 
turbulence (reference 35). 
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The previous conclusions depend upon the condition that the standard 
deviations of the amplitude modulated and mean value process are approxi- 
mately equal in atmospheric motions, that is, the value of the a 
parameter is approximately one. This condition is used as a preliminary 
estimate of the standard deviation of the mean value process. If measured 
atmospheric data show cases where the standard deviation of the mean value 
component is significantly different, then the previous conclusions must be 
revised accordingly. An associated question is the isotropy of the atmospheric 
motions. Examination of the measured data of references 18 and 19, which show 
the separation of the spectral functions into two fairly distinct parts, 
suggests that the mean value component process of the atmospheric velocity 
can be strongly anisotropic, giving significantly larger contributions to 
the horizontal velocity components. These data suggest that the a parameter 
can have significantly lower values for the lateral than for the vertical 
velocity component of atmospheric motion. The lateral response of the 
aircraft may accordingly be more significantly influenced by the mean value 
variation. 

The presence of the mean value variation in atmospheric motions has 
an indirect effect upon the calculation of aircraft response through the 
specification of the appropriate integral scale value for the local turbulence 
variations. The development of the integral scale values in measured 
atmospheric data requires the determination of two primary quantities: 
the integral scale of the mean value process and that of the local turbulence 
variations. Previous determinations of the integral scale value have not 
directly accounted for the mean value variations. This effect was included 
in the determination of the integral scale of the local turbulence variations, 
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that is, the local R component of the total process. This suggests that 
these determinations of the scale of turbulence, and. consequently the 
criteria value of 762 m (2500 ft), are too high. It is noted that 
suggestions for reducing the present scale of turbulence values have been 
made elsewhere (reference 35, for example). 

Development of Aircraft Design Criteria 

The application of the AMPM process to long-term measurements of 
atmospheric motions and the associated aircraft response requires the 
introduction of the concept of several types of turbulence. This concept 
was used in the original development of the amplitude modulated model 
(reference l) and has extensive experimental justification. The original 
random process is modified by introducing a conditional process, which 
is conditional on the type of turbulence. The probability density functions 
are weighted by the probability of the occurrence of each type. The quasi- 
steady exceedance expression of the modified AMPM process is obtained from 
the original expression, equation (A53). 
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Historically two types of turbulence are considered. These are often 
referred to as nonstorm (i=l) and storm (i=2) turbulence. There is an 
additional term, representing the probability of no turbulence, which is 
usually omitted in the literature. 

The exceedance expression of the modified process, equation (3^), 

can be applied to both the velocity components of atmospheric motion and 

the resulting dynamic system response under the quasi-steady assumption. 

For the case of an atmospheric velocity component, the deterministic mean 

value (w ) is zero and the two standard deviation factors (A and A ) 
o' ' r m 

have unit values by definition. 

The mean value variation introduces an additional atmospheric parameter. 
Consequently there are three atmospheric parameters for each type of 
turbulence: the probability parameter P , the amplitude intensity parameter 

b_^, and the mean value intensity parameter c_^. No data are available on 
representative values of the mean value intensity parameters since previous 
examainations of measured data have not considered the mean value variation. 
Also, these parameters cannot be obtained from long-term measurements of 
aircraft vertical acceleration, since the aircraft acceleration acts as a 
high-pass filter which largely removes the mean value variation. 

There remains the question of the effects of the mean value variation 
on the related problem of the specification of structural criteria for 
turbulence induced loads. The present approach is the specification of a 
maximum exceedance level for all aircraft loads. For the amplitude 
modulated process this requires the computation of two quantities: the 

standard deviation and the expected number of zero crossings of the 
response to the local R component process. From these quantities, the 
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deterministic mean value (one-g flight response) and a given set of 
atmospheric parameters, the exceedances of the response can he calculated. 
The introduction of the mean value variation requires the calculation of 
one additional quantity: the standard deviation of the response to the mean 

value process. This quantity can he determined hy the same computational 
procedures used for determining the response to the local K component 
process. Once the response quantities are known, the associated 
exceedances can he computed from equation (34) for a given set of 
atmospheric parameters. 

Aircraft design criteria are usually presented in terms of a gust 
intensity parameter. 
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For the amplitude modulated process the exceedance ratio depends solely 
upon U for given values of the atmospheric parameters. Consequently the 
aircraft design criteria can he specified as a required value of U^., 
which is a function, through the atmospheric parameters, of the altitude 
only. For the AMPM process this procedure is more complicated since the 
exceedance ratio of the response quantities depends additionally upon 
the a parameter. This is shown hy figures 7 and 8. Figure 7 shows a 
set of exceedance ratios as functions of the process level, and for 
several values of the a parameter. The exceedance ratios are computed 
for an altitude of 6.1 km (20,000 ft) using the atmospheric parameters 
of references 4 and 37 • No assumption on the values of the atmospheric 
parameters of the mean value process have been made, except that the a 
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parameters of the two types of turbulence are equal. From figure 7 the 
combinations of the gust intensity parameter U and the a parameter 
which correspond to a given value of the exceedance ratio are determined. 
These combinations are plotted in figure 8, which shows the Uq. and a 
parameter combinations which correspond to constant values of the exceedance 
ratio at the given altitude. It is noted that figure 8 also shows that 
the effect of the mean value variation upon the exceedance expression 
is very small if the a parameter of the response quantity has a value 
above two. 

The approach outlined above is based upon the exceedance ratio and 
not directly upon the exceedances , since the secondary effect of the 
expected number of zero crossings of the local E process for the 
response quantities is omitted. In order to use the number of exceedances 
directly, aJLl three parameters of the system response (A^, A^, N q ^) must 
be accounted for. The resulting exceedances can be computed from equation 
( 34 ). 



CONCLUSIONS 


The application of the AMPM ( amplitude-modulated-plus-mean) random 
process as an atmospheric turbulence model is examined and evaluated on 
a qualitative basis. The effect of the mean value variation is of primary 
interest since this is the distinction from the amplitude modulated process 
(the Press model, references 1 and 2) which is currently used to model 
atmospheric turbulence in aeronautical applications. It is concluded that 
the combination of amplitude modulation and a slow, random variation of the 
mean value is a better representation of measured properties of atmospheric 
motions. In particular, the AMPM process can account for the differences 
in the statistical properties of atmospheric velocity components and their 
gradients ; these differences cannot be accounted for by amplitude 
modulation alone. 

The correspondence between the properties of the AMPM process and the 
physical properties of atmospheric motions are examined. The structure 
of the random process suggests a possible correspondence with the structure 
of atmospheric motion: the lower frequency mean value process corresponding 

to the internal gravity waves (winds or drafts) and the higher frequency 
amplitude modulated process corresponding to the turbulence. The 
separation of the atmospheric motion into two elements appears in the power 
spectral density functions of atmospheric motions, which have the same 
general form as the functions of the AMPM process. 

The response of linear dynamic systems to the AMPM process is 
examined. Due to the structure of the process, the relative contribution 
of the mean value variation to the total system response is strongly 
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influenced by the dynamic properties of the response quantity. For 
response quantities which have the characteristics of a high-pass filter, 
such as aircraft accelerations, the relative contribution of the mean 
value variation is considerably reduced. For this reason, the gradients 
of the AMPM process show a larger deviation from Gaussian properties than 
the process itself. For response quantities which have significant 
response in the lower frequency range, such as aircraft displacements and 
velocities, the mean value variation can significantly increase the 
exceedances of the response at high response levels. Using 
estimates of aircraft dynamic properties, it is concluded that the 
effects of the mean value variation upon aircraft loads are small in 
most cases. However, the effects can be important in the measurement 
and interpretation of atmospheric motions. 



APPENDIX A 


SUMMARY OF MATHEMATICAL PROPERTIES OF AMPM PROCESS 

The mathematical properties of the amplitude-modulated-plus-mean or 
the AMPM random process are reviewed in this appendix. 

Process Definition and General Properties 

The AMPM process is formed from three independent random processes, 
which are identified as the local R, the amplitude S, and the mean 
value M components. The modulation of the local process hy the amplitude 
component forms the amplitude modulated process Z. 

z(t) = r(t) s(t) (Al)* 

The sum of the amplitude modulated and the mean value processes forms the 
AMPM process W. 


w(t) = z(t) + m(t) 


(A2) 


w(t) = r(t) s(t) + m(t) (A3) 

The properties of the AMPM process are determined hy the defining 
relation, equation (A3), and the properties of the three component 
processes, which are specified to he stationary and Gaussian with zero 
mean values. The notation for the variances of the component processes 
is 


*The numbering of the equations corresponds to that in reference 8; 
equation (Al) is equation (l) of that reference. 
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E[r 2 ] = A 2 

E[s 2 ] = b 2 (A10) 

E[m 2 ] = c 2 


The moments of the AMPM process are determined hy the defining relation; 
for example, the variance is 


_r 2 t .2 ,2 2 

E[w J = A b + c 


(A13) 


The ratio of the standard deviations of the independent amplitude modulated 
and mean value processes is a basic parameter of the AMPM process. 



(A14) 


The moments of the AMPM process can be expressed in terms of the a 
parameter; for example, the variance is 


E[v 2 ] = a 2 = e 2 (a 2 + 1) 


(A15) 


The fourth order flatness factor (or kurtosis) shows the dependence of 
the probabilistic structure of the AMPM process upon the relative 
contributions of the amplitude modulated and the mean value processes. 
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For a zero value of the a parameter, the flatness factor has the value 
of three for the Gaussian mean value process. For an infinite value of 
the a parameter, the flatness factor has the value of nine for the 
amplitude modulated process, the value of nine resulting from the product 
of two independent Gaussian processes. 

The autocovariance function of the AMPM process is related to the 
auto covariance functions of the amplitude modulated and the mean value 
processes. Using the defining relation and the independence property, 
the relationship is 

Y (t) = ¥ (t) + ¥ (t) (A22) 

w v ' z N ' nr ' 

where 

¥ (t) = T (t) 1' (t) 

X b 

The corresponding relation for the power spectral density function is 

§ (w) = $ (w) + $ (ai ) (A23) 

w z m 

The autocovariance functions can he expressed in terms of the (normalized) 
autocorrelation functions. 
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where 


V (t) = A 2 h 2 p (t) 

A A 

Y (t) = c 2 p (t) 

m K m 


(A29) 


P z (t) = p r (x)p s (x) 


Combining the previous relations and introducing the a parameter, the 
corresponding relations for the autocorrelation function and the normalized 
power spectral density function, which is the Fourier transformation 
of the autocorrelation function, are 
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(A31b) 


Quasi-Steady Approximation 

In application as an atmospheric turbulence model, the AMPM process 
is interpreted as the combination of a rapidly varying local component 
with slowly varying amplitude and mean value components. This concept 
leads to the quasi-steady approximation in which the dynamic properties 
of the amplitude process, and possibly the mean value process, are 
omitted in developing the dynamic properties of the AMPM process. The 
quasi-steady approximation presents a simple method for the analysis of 
the response of linear dynamic systems to the AMPM process. Assuming 
that the amplitude process affects the dynamic response of the system 
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only in a static manner, the dynamic response (or output) of a linear 
system to an AMPM process is also an AMPM process. 


v. (t) = r. (t) s. (t) + m. (t) 
m xn in in 


(A3' ) 


w (t) « r (t) s. (t) + m (t) 
out out v m out 
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Using the quasi-steady approximation, it is necessary to develop the 
response of the system to only the local R and the mean value components, 
which are stationary and Gaussian random processes. The variances of the 
three components of the input process are those of equation (A10) , except 
that a unit variance is specified for the input R component. The notation 
for the variances of the components of the output process is 


E[r 2 ] 
out 


E[m out ] 
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(A45) 
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The moments of the input and output processes are developed from equations 
(A3' ) and (a 43), and from the independence of the component processes. 


E[w 2 ] = c 2 (a 2 + 1) 

m in 


h[w 2 ] * A 2 c 2 (a 2 + 1) 
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The relation between the a parameters of the input and output processes 
is obtained from equations (Al4) and (a45). 


a + 
out 
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(a48) 


The values of the a parameters and consequently the probabilistic 
structures of the input and output AMPM processes can be different, 
depending upon the relative response of the d.ynamic system to the local 
E and to the mean value component processes. 

Another application of the quasi-steady approximation is the 
development of the derivative of the AMPM process. The exact derivative, 
which is developed from equation (A3), depends upon the derivatives of 
all three component processes. In the quasi-steady form, the derivatives 
of the amplitude and mean value components are omitted. The resulting 
quasi-steady derivative is • 


v (t) = r(t ) s(t) (A50) 

y. 

Quasi-Steady Exceedance Expression 

The exceedance expression (the expected frequency of the crossings 
of a given level) of a random process is developed from the joint 
probability density function of the process and its first derivative 
(references 38 and 39)- An analytical form can be developed for the 
exceedance expression of the AMPM process by omitting the derivatives 
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of the amplitude and mean value components. The resulting quasi-steady 
exceedance expression is 
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n A. 

W = — 

or 2tt A 

r 


The exceedance expression contains both the exponential dependence 
of the amplitude modulated process and the Gaussian dependence of the 
mean value process. The exponential dependence is dominant for large 
values of the a parameter. 

liffi N(v) = N e"l V l /Ab (A5T) 

cr* 30 or 

The Gaussian dependence is dominant for small values of both the ot 
parameter and the level of the process. 
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The exceedance expression approaches the exponential form of the amplitude 
modulated process for large values of the AMPM process. 


N(jw|»cct) « .N e' 


l/2a c 


I w I / ac 


(A59) 


Since the dynamic properties of the mean value component have "been eliminated 
by the quasi-steady approximation, the value of the expected number of 
crossings of the zero level becomes zero in the limit of small values of the 
a parameter, as indicated by equation (A58). The exact form of the 
exceedance expression, which includes the dynamic properties of all three 
component processes, is examined in appendix C. 



APPENDIX B. EFFECTS OF COREELATION BETWEEN THE AMPLITUDE 
AND MEAN VALUE COMPONENT PROCESSES 

In the formulation of the AMPM or total process, the amplitude and 
mean value component processes are specified to be independent. The 
resulting AMPM process is the sum of the amplitude modulated and the mean 
value processes which are independent. However turbulence theory 
introduces the concept of interaction between the two processes, specifically 
the flow of energy from the lower to the higher frequency process. This 
interaction raises the question of the effect of the assumed independence 
upon the properties of the resulting total process. In this appendix the 
formulation of a total process with correlated amplitude and mean value 
component processes is outlined. The effects of this correlation upon 
properties of the resulting total process are examined. 

The defining relation of the total process is 


w(t) 

= z(t) + m(t ) 

(Bla) 

z(t) 

- r(t) s ( t ) 

(Bib) 


The amplitude and the mean value components are specified to be 
correlated Gaussian processes. 

E[sm] = p be (B2) 

sm 
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The local R component is specified to he independent of the other two 
component processes . The resulting total process is the sum of the 
amplitude modulated and the mean value processes. The relation between 
these two processes is shown by their joint moments, for example, 


E[zm] = 0 

(B3a) 

t) r 2 2-i . 2 , 2 2/_ n 2 ■. 

ELz m ] = A b c (1 + 2p ) 

(B3b) 


Equation (B3a), which follows from the independence and zero-mean 
properties of the R component, shows that the amplitude modulated and 
the mean value processes are uncorrelated. Equation (B3b) shows that 
the amplitude modulated and the mean value processes are not independent 
(for a non-zero value of the correlation coefficient). 

The moments of the corresponding total process are developed from 
the defining relation, equations (Ela) and (Bib). All odd order moments 
are zero. The relations for some of the even order moments are: 

E[w 2 ] = c 2 (a 2 + 1) (BUa) 

E[w^] = 3c^[3cx 4 + 2a 2 (1 + 2p 2 ) + 1] 

sm 

E[v 6 ] = 15c 6 [15a 6 + (9a^ + 3a 2 )(l + 4p 2 ) + l] 

sm 


(B4b) 
(b4c ) 
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Examination of the basic form of the moments shows that the correlation 
coefficient always appears as an even power. The correlation is thus a 
second order effect for small correlation values. 

The effect of the correlation is also shown by the associated fourth 
order flatness factor. 

3[3a^ + 2a 2 (l + 2p 2 ) + 1] 

M V(W) a (B5) . 

(a + 1) 

The correlation increases the value of the flatness factor, thus showing 
tendency away from a Gaussian distribution toward that of the amplitude 
modulated process. However, the minimum value of three and the maximum 
value of nine for the flatness factor are not changed by the correlation. 
In order to estimate the numerical effects of the correlation upon the 
flatness factor, the case of equal values of the standard deviations of 
the amplitude modulated and mean value process is considered. 

M>(wi a = 1) * 4.5 + 3 p 2 (B 6) 

H SDl 

The correlation coefficient must have an appreciable value in order to 
significantly change the value of the fourth order flatness factor. 
Conversely, if the amplitude and mean value component processes are 
weakly correlated, it will be difficult to determine the extent of 
correlation from measured data by use of the flatness factor. 
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The relations for the characteristic and probability density functions 
of the total process follow from the defining relations, equations (Bla) 

» and (Bib). The relation for the characteristic function is 

C (6) * /“ f° e imQ C (s0) p (s,m) ds dm (b?) 

Since the required functions of the three Gaussian component processes 
are known, the characteristic function of the total process can be 
determined. 

The effect of the correlation upon the exceedances of the total process 
is an important question for aeronautical applications. The exceedance 
expression is developed from the joint distribution of the process and its 
first derivative. Using the quasi-steady approximation (reference 8), 

B(v) = f m N r (~ 7 ~~) P(s,m) ds dm (B8) 

Equation (B8) is integrated numerically since the required integrations 
appear to be intractable. Figure 9 shows the resulting exceedance 
expression as a function of the ratio of the level of the total, process 
to the standard deviation (Ah) of the amplitude modulated process. The 
exceedance curve is shown for an a parameter of . 60 , which is 
approximately the value for the maximum differences \ the effects of the 
correlation upon the exceedance expression vanish in the limits of zero 
and infinite a values. The correlation has the effect of increasing 
the exceedances at high values of the process level. For a value of the 
non-dimensional process level equal to six, which is about the critical 


level in aeronautical applications , the exceedances are increased by 
less than 20% for a correlation coefficient of .30. 

In summary, the introduction of correlation between the amplitude 
and mean value component px-ocesses eliminates the independence property 
of the amplitude modulated and mean value processes. Thus the two 
processes can interact as required by basic turbulence theory. If the 
amplitude and mean value processes are weakly correlated, then the 
correlation is a second order effect; the effects of the correlation 
on the fourth order flatness factor and on the exceedance expression 
of the total process are fairly small. 



APPENDIX C. COMPARISON OF THE EXACT AND QUASI-STEADY 
EXCEEDANCE EXPRESSIONS 


The exact exceedance expression of the AMPM or total process is 
examined and compared -with the quasi-steady form of the exceedance 
expression, equation (A54), in this appendix. In the quasi-steady 
approximation the dynamic properties of the local R component process 
are considered, -while those of the amplitude and mean value component 
processes are omitted. In the exact form of the exceedance expression 
the dynamic properties of all three component processes are considered. 
The present development extends that of references 11 and 12. Although 
the exceedance expression is examined in general form, the modulated 
form of the total process is of primary interest, that is, the amplitude 
and mean value component processes are slowly varying relative to the 
local R process, with the quasi-steady approximation being a limiting 
case. 

The exceedance expression is developed from the joint probability 
density function of the random process and its first derivative 
(references 30 and 39)> 

N(v) = / w p(w, w) dw (Cl) 

The indicated integration operation essentially gives one-half of the 
first absolute moment of the first derivative. (This moment is assumed 
to exist in the subsequent development.) For the total process it is 
convenient to develop the exceedance expression from the joint character- 
istic function of the process and its first derivative. 


53 



Tne absolute moments of a random variable are developed from the charac- 
teristic function (reference 40). 


E[ | x j n ] = /“ |x[ n p(x) dx 


1 r 00 3 

_ . n+1 o „ >.n 
2m 3? 


tc(5) + c (-?).] 


d5 


(C3) 


where n = odd. 

The relation for the Fourier transformation of the exceedance 
expression is obtained by combining the previous relations. 


F{N(w) } = 


-1 3 

7T O 30^ 


[ c (0 v> 0 d )] 



(c4) 


wnere C( 0 , -0 ) = C( 0 0 ) 

' w d' ' w d 

Introducing the inverse Fourier transformation, the relation for the 
exceedance expression is 


= v Ol sr f c < 9 

7T d 


w* 6 d^ 


d6 

cos(w0 )d0 — r— 

w W 0J 


(C5) 


The joint characteristic function of the total process and its first 
derivative is 



55 


obtained from equation (C?)» Since the total process is the sum of the 
independent product and mean value processes , the characteristic function 
is the product of the functions of these two processes. 


C .(0 , 0 ) = C .(0 ,8 )C .(0 , 0 ) 
wv d zz w d mm w d' 


(c 6) 


The joint characteristic function of the product process Z and its 
first derivative is given in reference 7* Since the three components 
of the total process are stationary and Gaussian, it is convenient to 
express the standard deviations of their first derivatives in terms of 
the expected number of positive slope zero crossings. The appropriate 
notation for the local R component process is 


N 


i A. 
_ _jL _r_ 

or 2 tt A 


(C7) 


where 


ft 2 -nr 2-, 

A r = E[r ] 


A? = E[f 2 ] 


The notation for the amplitude component process is 


ST 


os 


r • 2 -i A 
ELs~~] r 


£ s U Qr \|E[s2] A f 


(C8) 


The notation for the mean value component process is 
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e - 


om 


!l£j A r 


i tr r 2 i a. 

or 1/ E[m j . r 


<C9) 
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The concepts of amplitude modulation and the slow variation of the mean 
value require that the e parameters have values much less than one, 
with zero values corresponding to the quasi-steady approximation. 

The joint characteristic function of the total process and its first 
derivative is obtained hy combining the previous relations. 


c ( 0 ,e,) = [u 2 x 2 + 

w s a' u 


/ 2 2w n , 2 2 2^i 

(1 + a y ) (1 + e s a y ) ] 


- 1/2 


• exp(- x 2 



2 2 , 

e y ) 

m 


(CIO) 


where 


x = A b0 /a 
r w 

y = A r l30 d /oi 


The relation for the exceedance expression is obtained by combining the 
joint characteristic function and equation (C5). 


N(w) 2 „ n -l/2 , 1 2 1 2 2 v 

= ~ fl> exp(- ^ x “ 2 e - y ) 
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f = e 2 + [1 + e 2 (1 + 2a 2 )y 2 ]h _1 

m s 


(Oil) 
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In the special case of both £ parameters equal to zero the integrals 
reduce to a standard form., giving the quasi-steady exceedance expression, 
equation (A5H). In the general case the integration appears to be intractable 
and is evaluated numerically. 

The exact exceedance expression is plotted in figures 10 and 11 as a 
function of the level of the total process (ratioed by the standard 
deviation, Ab, of the product process), and for several values of the 
parameter. The dynamic properties of the amplitude component process have 
been eliminated by setting the value of the £ g parameter to zero. Figure 
10 shows the exact exceedance expression for the value of the a parameter 
equal to one, in which case the variances of the product and mean value 
processes are equal.. Figure 11 shows the exact exceedance expression for 
the value of the a parameter equal to .5, in which case the mean value 
component is the dominant contributor to the total process. The exceedance 
expressions for the quasi-steady case correspond to the zero value of the 
£ parameter. As the value of that parameter is increased, the values of 
the exceedance expression are increased at all levels. Thus the dynamic 
properties of the mean value process uniformly increase the exceedances 
of the total process. The relative increase is greater for smaller values 
of the a parameter, that is, for larger static contributions of the mean 
value component to the total process. The largest relative increase in 
the exceedances occurs at the zero level of the total process. The exact 
and quasi-steady exceedance expressions are equal in the limit of large 
values of the level of the total process. The results of figures 10 and 11 
show the following relation for the special case of equal values of the 
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expected number of zero crossings of the local R and mean value component 
processes, that is, the £ m parameter equal to one. 


H(w = 0; a, E = 1, £ 
m s 


0) * N 


or 


N 

om 


(C12) 


This relation follows from the general expression, equation (Cll), which 
can he integrated analytically in this special case. 

The concept of the AMPM process specifies that the mean value process 
is slowly varying relative to the local R process. Consequently the 
value of the £^ parameter is generally much less than one. However there 
can he exceptions to this which are within the concept of the AMPM process. 
One example is the response of a narrow hand-pass system to the AMPM 
process. The frequency content of the system response, including hoth the 

amplitude modulated and mean value processes, will he dominated hy the 

frequency of the narrow hand-pass system. Consequently the expected 
number of zero crossings of the local R and the mean value processes of 

the response will he almost equal; the value of E^ parameter of the 

system response process will he approximately one. 

The effects of the dynamic properties of the amplitude process, which 

are expressed through the £ g parameter, equation (C8), were examined hy 

numerical integration of the exact exceedance expression. Numerical 

results show the following general effects for values of the £ g parameter 

less than .25. If the value of the a parameter is one or greater, the 

effect of the £ parameter upon the exceedances is always less than 
s 
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The effect is larger for smaller values of the a parameter. The effect 
of the e parameter can significantly increase the exceedances in the 
of very small values of both the a anu the parameters. 


case 
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Figure 1. Example of atmospheric velocity components showing random amplitude and 
mean value variations (flight 32, run 2 from NASA-Langley MAT program, reference 14) 




log frequency 


Figure 2. General form of the two terms of the power spectral density 
function of the AMPM process 
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Figure 3- Exceedance .ratio as a function of the process 
level non-dimens ionalized for 4 the AMPM process, and for 
several values of the a parameter 
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Figure 4. Exceedance ratio us a function of the process 
level non-dimensional ized for the amplitude modulated process 
and for several values of the a parameter 
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Figure 5 . Exceedance ratios for the quantities of a first- 
order differential equation with the parameters of equation (29) 
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Figure 6. General relation between the power spectral density functions 
of the two terms of the AMPM process 
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'Figure 7. Exceedance ratio as a function of process level, 
for several values of the a parameter (6.1 km altitude) 
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